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We analyse dark energy models where self-interacting three-forms or phantom fields drive the 
accelerated expansion of the Universe. The dynamics of such models is often studied by rewriting 
the cosmological field equations in the form of a system of autonomous differential equations, or 
simply a dynamical system. Properties of these systems are usually studied via linear stability 
theory. In situations where this method fails, for instance due to the presence of zero eigenvalues in 
the Jacobian, centre manifold theory can be applied. We present a concise introduction and show 
explicitly how to use this theory in two concrete examples. 
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; I. INTRODUCTION 

O ! Observational cosmology provides strong evidence that the Universe is currently undergoing an epoch of accelerated 
expansion (see e.g. The driving force of this expansion is called dark energy, and it appears that the Universe's 
■ energy content is made up of 74% of this mysterious dark energy. The second biggest energy content of the Universe 
is a non-baryonic matter, interacting only via its gravitational attraction, called dark matter. The existence of such 
a substance has long been known. It is implied by the flattened galactic rotation curves observed by Zwicky as early 
,; as 1933. Neither dark energy nor dark matter have been detected directly. 
^ ' The simplest model of dark energy is the cosmological constant A, originally introduced by Einstein in 1917 in 
C — ■ ', order to construct the Einstein static universe but later abandoned. The cosmological constant appears to be in very 
\l good agreement with observational data, however, its physical interpretation is unsatisfying when adopting a particle 
physics point of view. In this context it is interpreted as a measure of the vacuum energy density which leads to the 
well known cosmological constant problem 0]. Another important issue is the question why the energy contents of 
dark energy and dark matter are of similar order of magnitude today. This requires finely tuned initial conditions 
from which the Universe must have started its evolution. 

Some of these issues can be addressed by introducing a dynamical model for the evolution of the dark energy. For 
instance, one could model dark energy as an evolving scalar field with an asymptotically flat potential. If dark energy 
interacts with dark matter, the cosmological coincidence problem may be alleviated. Therefore, various interacting 
dark energy models have been proposed and investigated [3142 Ol| . Some models such as quintessence have been studied 
in It is logical to also consider more complicated fields, like spinors 122 
' fields. One such suggestion was recently studied by Koivisto and Nunes [2/ 

studied their cosmology. Those fields can yield a viable cosmology where they drive infiation in the early time and 
become the driving forces of the accelerated expansion of the Universe at late times. This has also been investigated 
in 

Besides the exact form of the field describing the dark energy, one can exploit another freedom in current cosmo- 
logical models. Since both dark energy and dark matter are not understood fundamentally, there are no a priori 
constraints on interactions between these components. In all these scenarios it seems relevant to consider some degree 
of interaction (other than the gravitational one) between dark matter and dark energy. Dark matter is primarily 
needed in the early Universe to support the formation of structure. On the other hand, dark energy is required to 
drive the late time accelerated expansion of the Universe. Thus, it seems plausible to study models with an energy 
transfer from dark matter to dark energy as time progresses. This idea has motivated a vast amount of literature 
on this topic, see for instance [l^l and references therein. Since there are no selection criteria leading to a specific 
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coupling, any coupling that might be studied will necessarily be phenomenological, with some models having a good 
physical justification while others are motivated by mathematical simplicity. 

II. DARK ENERGY MODELS 
A. Basics 

We work in a flat FLRW spacetime whose line element is given by 

ds^ = -dt^ + a^{t){dx^ + dy^ + dz^). (1) 

Neglecting radiation and baryons, the conservation equations for the background dark matter (DM) fluid and the 
dark energy (DE) are 

PDM = -iHp-DM + Q, (2) 
/ODE = -3i7(l + wde)/3de - Q, (3) 

where Q is an arbitrary coupling and the subscript DE stands for a generic dark energy model to be specified. The 
conservation equations are subject to the Friedmann constraint 

^2 

= y(/5DE + PDM), (4) 

where = SttG/c^ is the gravitational coupling strength. 

B. Three-forms 

Let us now consider the case where a three-form field generates both early (infiation) and late time (dark energy) 
acceleration. The non-zero components of the most general three-form field in this geometry is described by j3, [111 

A,jk^a{tfe,juX{t), (5) 

where X{t) is a scalar function of time and Cyfe is the standard three-dimensional permutation symbol. The equation 
of motion of this field is given by a modified Klein-Gordon equation, which reads 

X = -iHX ~Vx- 3HX - -. — . (6) 

x + 'mx 

Following [13, HI] , we choose an exponential potential 

F(X) = Vbe-^^, (7) 
where A is a dimensionless parameter and Vq > 0. Using this potential, the field equations become 

Px^l{X + 3HXf + V{X)-Q, (8) 
Px^-^{X + 3HXf-V{X) + VxX, (9) 

= ^{^iX + 3HXf + ViX)+pDM), (10) 

H = ^{VxX + py,m)- (11) 

To construct a dynamical system, we use the following dimensionless variables 

1 



V6H 



{X + 3HX), (12) 

(13) 



— arctan^^, (14 
TT ^/6 

PDM 



(15) 
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Note that we have now deviated from [2j] . Our choice of variables has the advantage that the phase space is compact 
by construction and we do not have to worry about the presence of critical points at infinity. These variables are 
motivated by noting that the Friedmann constraint now becomes 



+ + = 1. 



(16) 



Note that by construction — 1 < z < 1, and moreover — l<a;<l,0<y<l and < s < 1. Our phase space is 
therefore a half cylinder of height 2. 

The equation of state parameter for the three-form field is defined by wx = Px / px and thus can be written as 



wx 



-1 



VxX 



-1 



Px + y 

Similarly, the total equation of state parameter becomes 

wtot = -x'^ - 2 2/^ fs + \/6Atan 



nz 

Y 



TTZ 

T 



(17) 



(18) 



Recall that the condition for acceleration is wtot < —1/3. We choose the coupling Q to be of the form Q = apcH . 
This results in the following autonomous system of differential equations 
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(19) 

(20) 
(21) 



Note that the equation system is invariant under the map y ^ —y and thus we can constrain the analysis to the 
y > case. This is in fact not surprising since our potential is positive definite. The number of critical points of this 
dynamical system depends on the coupling parameter a. Starting with a = 0, we arrive at the results of [2^ which 
are 
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TABLE I. Critical points of the uncoupled model. 



In the presence of a coupling a 7^ 0, the number of critical points changes. Specifically, the point B splits into 2 
different critical points, see Figure [TJ The points A± remain unchanged, however, their eigenvalues do change. This 
is summarised in Table HIl 

As the coupling strength increases to its maximally allowed value, a — > 3, the two points B± move towards the 
points A±. When a = 3 these points merge and the system has two critical points, each with two zero eigenvalues. 
Note that we do not analyse this degenerate case. 



C. Phantom dark energy 

Phantom dark energy models have been of interest since their first introduction by Caldwell [HI. Cosmological 
observations place the dark energy equation of state close to —1, however, this value could also be approached from 
below, a case which is not excluded by observations. In a recent paper, Leon & Saridakis [3C| have considered a 
varying-mass model for dark matter particles in the framework of phantom cosmologies. They considered a phantom 
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TABLE II. Critical points of the coupled model. 



dark energy model with power-law potential interacting with dark matter. Dark energy is modelled as a scalar field 
with negative kinetic energy whose energy density and pressure are given by 

P4.^-\\? + V{<i>), (22) 

V4>^-\^^ ~y{^)- (23) 

We denote the equation of state parameter by = p^j p^. Note that the sign of the phantom kinetic term is opposite 
to that of an ordinary scalar field. The dark matter energy density is assumed to depend on the mass of the dark 
matter particle which in turn is assumed to depend on the field pdm — -^^dm(0)"-dm where tidm is the number 
density which is determined by jt-dm + 3i?nDM — 0. The field dependence of the mass can be interpreted as a coupling 
between the dark energy and the dark matter. Therefore, the dark matter energy density satisfies the evolution 
equation Q where the coupling Q is determined by the dependence of the mass on the field 0. This means 

din Mdm(0) 

Q = 77 <P/3DM- (24) 

d(f) 

Assuming that the total energy-momentum tensor is conserved, the dark energy will satisfy equation ([3]). Both 
evolution equations will be subject to the Friedmann constraint Q. 
Similar to P^ -([T5 |) . we define dimensionless variables as follows 

(25) 



Veil' 



It is useful to express the density parameter and the equation of state in form of these variables, which gives 

^,-^ = --' + y', (28) 

9 2 

— X — y 

^'t' = — — 2^ (29) 
wtot = -x" - y^. (30) 
In these variables, the cosmological field equations take the form of the following dynamical system 

x' = -3x + \x{l -x^-y^)-^l- ^z(l + x^~ y% (31) 

v'-\v{^-x'~v')~^, (32) 
z' = -xz^. (33) 
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TABLE III. Critical points of the phantom dark energy model. 



This system possesses two physically meaningful critical points which are non-hyperbolic since there exists at least 
one zero eigenvalue in each of them, see Table IIIII 

Thus, centre manifold theory is required in order to study the nature of these critical points as remarked in [30| . 
Note that point A is always unstable since we have one positive and one negative eigenvalue. 



III. CENTRE MANIFOLD THEORY 

A. Mathematical background 

When discussing the mathematical background of the centre manifold theory, we closely follow Carr [3l| and 
Wiggins fs^]. In the presence of zero eigenvalues in the stability matrix at the critical point, linear theory fails to 
provide information on the stability of that point. The main aim of the centre manifold is to reduce the dimensionality 
of the system near that point so that stability of the reduced system can be investigated. There always exists an 
invariant local centre manifold passing through the fixed point to which the system could be restricted to study its 
behaviour in the neighborhood of the fixed point. The (in)stability of the reduced system determines the (in)stability 
of the system at that point. 

Let a; € M*^ and y £ . An arbitrary dynamical system with zero eigenvalues in the stability matrix can always be 
written in the following form 

Ax + f{x,y), 

y = By + g{x,y), (34) 

where 

/(0,0) = 0, Df {0,0)^0, 

5(0,0) = 0, i^g(0,0)=0. (35) 

Here we assume that the critical point is located at the origin and Df denotes the matrix of first derivatives of the 
vector valued function /. A is a c x c matrix having eigenvalues with zero real parts and _B is an s x s matrix having 
eigenvalues with negative real parts. 

Note that a dynamical system with zero eigenvalues can always be rewritten into the above form by virtue of a 
linear change of coordinates. We will show this construction explicitly when we discuss the three-form and phantom 
models below. 

Definition. We call the space 

W'iO) = {{x,y) G R= X W\y = h{x),\x\ < S,h{0) = 0,Dh{0) = 0}, (36) 
for 6 sufficiently small, the centre manifold for the system (l34l) . 

One can think of this manifold as the space where the stable directions y are parametrised by the unstable directions 
X. One can now consider the dynamics of the system restricted to this manifold [3l[ . Since the y = h{x), one is left 
with the reduced equation 

■ii^ Au + f{u,h{u)), weR^ (37) 

for sufficiently small u. The (in)stability properties of the reduced system (1571) imply (in)stability properties of the 
fuU system ((M)) . 
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Next, we need to construct this centre manifold explicitly. By differentiating the defining equation y = h{x) with 
respect to the independent variable we get y — Dh{x)x where we used the chain rule. Eliminating x and y via (j34p . 
one arrives and the following quasilinear partial differential equation which h has to satisfy 

M{h{x)) = Dh{x) [Ax + f{x,h [x))] ~ Bh{x) - g{x, h{x)) = 0, (38) 

from which h can be found and inserted into (j37p to study the reduced system. In general it will not be possible to 
solve (|38p analytically. However, such knowledge is in fact not needed since we are only interested in the reduced 
system near the critical point. Therefore, it suffices to simply assume h to be of the form h[x) = ax^ + bx'^ + 0{x'^) 
and to determine to first few non-trivial terms of the Taylor expansion of h. Note that this is always possible and will 
result in a unique solution. 



B. A simple example 

Let us start with a simple example [32.] that shows the above construction. Let {x, y) G and consider the system 

x = x^y-x^, (39) 
y = -2/ + a;2, (40) 

whose only critical point is {x,y) — (0,0) with eigenvalues and —I. 
Therefore, in this example one has 

A = Q, B = -l, (41) 

f{x,y) =x'^y-x^, g{x,y)^x^. (42) 

Assuming that h{x) is of the form 

h{x)^ax^ + bx^ +0{x'^), (43) 

we find that 

Af{h{x)) = {2ax + 3bx^ + 0{x^)){x^ {ax^ + bx^ + 0{x'^)) - x^) + {ax^ + bx^ + 0{x^)) - x^ (44) 
^ {a - l)x^ + bx^ + 0{x^) ^ 0. (45) 

Comparing coefficients, we find a = 1 and b = and thus the centre manifold is given by h{x) ~ x'^ + 0{x^). 
Now, the evolution equation restricted to the centre manifold (I57| is given by 

x^x^ + Oix"^), (46) 

whose dynamics are such that for x sufficiently small, x = is unstable. 

We would also like to note that it is not allowed to simply approximate the centre manifold by y = 0. Putting 
this into the equations (I40p . we would arrive at the reduced equation x = —x^ which would indicate stability. More 
details can be found in j31j, i32J 



IV. APPLICATIONS TO DARK ENERGY MODELS 



A. Three-forms 



Applications of centre manifold theory to cosmological models have been studied previously, see e.g. [33|,|3J]. We 
now apply this to the system (|19p -(|2T |) for the critical point In order to do so, we must transform this system 
into the form of (|34p . Firstly, we introduce new coordinates X~x — 1, Y~y and Z = z — 1/2 which move the 
point (1, 0, 1/2) to the origin of the phase space. By computing the matrix of eigenvectors of the stability matrix of 
the system in X, Y, Z, we introduce another set of new coordinates 

1 o\ /x\ 

-3/(™) 1 y\. (47) 
3/(7ra) o) \Z ) 
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In these coordinates, our system of equations is now in the correct form 



B 



-3 



/ A 

u 


J 


V 




V"/ 




3rm 




, we 


find 







-3 + 


a / ' 




fu] 




^ non 


V 




hnear 


[w) 




^terms 



1 3 1 

/ = 7ra(6 + V&X)uw lY^a^uiD^ 

6 2 6 

_ ( {a + 2a27ra - 3)^2 + Sagu^ + 0{u^) \ 
^ ~ I (3 + 2&2vra)u2 + 63(3 _ a)u^ + ©(w") ) 



-uA ( 1 + (u- 



1) tan 
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(48) 
[v, w) is a 

(49) 
(50) 
(51) 



Note that knowledge of g up to order 0{u^) is sufficient to construct the centre manifold up to the desired order. The 
complete expressions for the components of g can be found in Appendix |^ 
The centre manifold can now be assumed to be of the form 



a2U^ + 03"^ + 0{u^) 



It has to satisfy equation 



which explicitly reads 



1 



3{a + 2a2TTa - 3)u^ + 3a3U^ + 0{u'^), 



2Tia \^(3 - a)(3 + 263^)^2 + 63(3 - a)u^ + 0{u^),^ 

Solving for the four constants 02,03,62 and 63, we obtain 

3 — a 



= 0. 



02 = 
h2 - 



2'Ka ' 
_-3_ 

27rck ' 



03 = 0, 
63 = 0. 



We can now study the dynamics of the reduced equation (|37|). which becomes 



(52) 

(53) 

(54) 
(55) 

(56) 



Therefore, we find that the point is stable according to centre manifold theory. When we repeat this calculation 
for A_ we find the opposite results, namely, this point is unstable. Figure [T] shows the phase spaces of the coupled and 
the uncoupled models. Note that this information is not obvious from the trajectories in the phase space diagrams 
since one is tempted to conclude that both and A_ are stable. 



B. Phantom dark energy 

We extend our work to a scenario of varying-mass dark matter particles in the framework of phantom cosmology, 
and our results obtained by performing centre manifold analysis are consistent with [j30|. 

We firstly look at model 2 of jSOj ■ The autonomous system of differential equations is given by 

x' = -3x + ^x{l -x^-y^)-^- ^z{l + x^~ y^), (57) 

y 2^*- ~y > — 2~' ^ ^ 

z' = -xz^. (59) 
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The critical points are (0, 0, 0) and (0, 1, 0), which are both non-hyperbohc. However, the former has a positive and 
a negative eigenvalue, from which it can be deduced the point is unstable. As for the second point, there are two 
negative eigenvalues and a zero eigenvalue. Therefore, linear theory fails to provide information about the stability 
of that point. It is required to perform centre manifold analysis in order to study the nature of the stability of the 
second point. We first transform the coordinates into a new system X = x,Y = y— 1 and Z — z, such that the 
critical point is at the origin. Next, we introduce yet another new coordinate system 



(60) 





(o 
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fx 
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f) 




[wj 


lo 










In these new coordinates the equations are transformed into 



As before, we compare this with the general form 
a two-vector component vector. Thus, we obtain 




(non \ 
linear 
terms / 



(61) 



and deduce that x = u is a scalar function while y = (v, w) is 



A = 0, B = 




f 
9 



;U^{6v — uX), 



(62) 
(63) 
(64) 
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where we give the exphcit expressions of the components of g in Appendix |B] 
The centre manifold can now be assumed to be of the form 

[b2U^ + b3U^ + 0{u'^) / ■ 



Equation (j38|) becomes 



where the 02 , 03 , 62 and 63 are computed as 



(65) 



fSa2U^ + jk4^^432a3 + 72b2iX-2^i)-X\A + X-2^i))+Oiu^)\ 



02 = 0, 03 = — , (67) 

62 = ^, 63 = 0. (68) 

The dynamics of the system restricted to the centre manifold is given by 

u=-u^ + 0{u^). (69) 
6 

Thus it is clear from this equation that the point is stable if A < and unstable if A > 0. This result is consistent 
with [30] where the same result was obtained by using the method of normal forms. 



V. CONCLUSIONS 



One of the main aims of this paper was to give a concise introduction into centre manifold theory, which is useful for 
studying dynamical systems with zero eigenvalues in the stability matrix. While this method has recently been applied 
to some models, a simple introduction with applications to cosmology does not seem to appear in the literature. We 
provide such an introduction and point out the necessity of using centre manifold theory. 

We then apply this method to two dark energy models, one based on three-forms and the other one based on 
phantom dark energy. For the phantom dark energy model we confirm previous results and show their derivation 
explicitly. For the three-form model we can conclude that only one of two potentially stable points is in fact stable, 
a result not accessible from linear stability theory. Moreover, the trajectories of the phase space seem to suggest 
stability of the unstable point, emphasising that critical points with zero eigenvalues need to be investigated carefully. 

In order to study physically relevant dynamical systems derived from cosmological field equations, it is important 
that we do not restrict our models based on simplicity and the ability to use linear stability theory. We should 
aim towards constructing the most physically motivated models and study them irrespective of the mathematical 
challenges involved. There exists a huge amount of mathematical literature on dynamical systems, much of which 
still awaits to be applied to physics. 
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Appendix A: Explicit formula of the three-form model 

The two components of g in the three-form model are given by 



1 /9 27 6 

gi — - [ — h 18(w + w) + 3wa(l -I- Sttw) + ir^w'^a'^ H — (3 cos[7r(?j -I- w)] — (3 -I- ■kwo) sin[7r(i; + w)]) 

6 \7r 7r(3 -|- TTwa) tt 



1 , / 27 3 

6 \ 7r(3 -I- TTwa) ira 



3 (3 - Vga) + \/6(3 + 



1 /I 

— TT \- V + W 

2 V2 



, (Al) 



92 



3u2 ((3 - QTTw)a - n'^w'^a^ - 9) 



27rQ:(3 + ■Kwa) 



— \ —u^X ( 3 — (3 + TTwa) tan 
na \ 2 ' 



— TT 1 V + W 

2 \2 



TTw'^a (27 + 3(1 + 4:TTw)a + n'^w^a'^) 
6(3 -I- TTwa) 



(A2) 



Appendix B: Explicit formula of the phantom model 

The two components of g are given by 

3 1 

= --V (z;2 + w{2 + w)) + -u {w{2 + w)[2^ - A) + w^(3A - 2/x)) 



^^u\X{A + 3A - 4m) + Y^"'A^(4 + A - 2/i), (Bl) 
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